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ANALYSIS OF OSO DATA TO DETERMINE THE STRUCTURE 


AND ENERGY BALANCE OF THE SOLAR CHROMOSPHERE 
Grant NGR 09-015-210 

The objectives of the research supported by this grant are twofold: 

1) The determination of an improved empirical model of the temperature-density 
structure of the solar atmosphere based on OSO 4 and 6 data* and 

2) The determination of an associated theoretical model based on radiative 
balance, thermal conduction, and shock heating. 

The following work has been performed to accomplish these objectives. 

We have carried out a detailed reexamination of the temperature-density structure 
of the photosphere and low chromosphere. The middle and upper chromosphere, which 
directly emits most of the OSO spectrum, is sensitive to conditions in this underlying 
region of the atmosphere. Our model of this region is based on a unified compilation 
of all recently published broadband flux and central intensity observations of the solar 
spectrum from 500 pm in the far infrared to 1220 A in the far ultraviolet. This exten- 
sive compilation includes the OSO 4 and 6 observations in the wavelength range 1400 to 
1220 A. A complete report on this work by Vernazza, Avrett, and Loeser has been 
widely distributed in preprint form and will appear in The Astrophysical Journal 
Supplement Series . 

We have also calculated a new model of the quiet solar atmosphere in the height 
range between the temperature minimum and the upper part of the chromosphere- 
corona transition region. This model is based on statistical equilibrium calculations 
of H, He I, He II, Si I, Cl, and other ions. The results of these calculations will 
appear in forthcoming papers to be submitted to The Astrophysical Journal. 
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Finally, we have studied the propagation and dissipation of acoustic waves in the 
solar atmosphere and the dependence of the dissipation on the energy flux and the 
period of the waves. We conclude that the low chromosphere is heated by short-period 
waves. An account of this work was presented by Kalkofen at a recent solar physics 
meeting of the American Astronomical Society. The text of that paper follows as 
Appendix A. 

In addition, we reproduce in Appendix B the text of the paper "Formation of the 
Solar EUV Spectrum, " presented by Avrett at the same meeting. 
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HEATING OF THE LOW SOLAR CHROMOSPHERE 


Wolfgang Kalkofen and Peter Ulmschneider 
INTRODUCTION 

The temperature of the solar atmosphere decreases monotonically from the 
interior to the height of the temperature minimum and then rises again into the chromo- 
sphere and corona. This temperature rise is due to nonradiative energy input into the 
chromosphere. Biermann (1948) and Schwarzschild (1948) have suggested that the 
energy that heats the chromosphere is carried in mechanical waves that are generated 
in the convection zone and dissipate their energy in the chromosphere. Osterbrock 
(1961) has estimated the flux, and Stein (1968) has calculated frequency spectra of the 
acoustic waves. In this paper, we want to study the propagation and dissipation of 
acoustic waves in the solar atmosphere and the dependence of the dissipation on the 
energy flux and on the period of the waves. 

Observations place two constraints on the energy flux and the period of the waves: 
They must be such that the height at which the sound waves steepen into shocks and 
dissipate a large fraction of their energy agrees with the location of the temperature 
minimum and that the energy flux remaining in the waves thereafter is sufficient to 
account for the radiative losses of the chromosphere. 

The theory of wave generation gives the energy flux and the spectrum. These 

calculated quantities are very uncertain because of the dependence of the flux on the 

eighth power of the velocities in the convection zone, which are only poorly determined 

from mixing-length theory. Osterbrock (1961) lias determined a value for the noise 
7 2 

flux of 3 X 10 ergs/cm sec; Stein's (1968) value is slightly larger. For several 
turbulence spectra, Stein has calculated the corresponding noise spectra, which show 
broad maxima with periods between 20 and 80 sec. 
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The dissipation of sound waves in the low chromosphere has been studied by 
Ulmschneider (1970). He has concluded that the low chromosphere must be heated by 
short-period waves whose periods are on the order of 30 sec. This conclusion is sup- 
ported by Stein's spectrum of the acoustic flux. It is not supported by observations, 
but there are also no observations that would rule out heating by short-period waves. 

7 

Acoustic waves with periods near 30 sec and an initial energy flux of 3 X 10 ergs/ 

2 

cm sec at the bottom of the photosphere dissipate their energy just beyond the tem- 
perature minimum region. Waves with longer periods, and therefore longer wave- 
lengths, travel farther into the chromosphere before they dissipate, and waves with 
a shorter period dissipate earlier. If waves with periods much below 30 sec carried 
enough energy, they would cause the chromospheric temperature rise at lower depth; 
and thus the temperature minimum would lie at a lower depth, contrary to observation. 

RESULTS 

We have solved the time -dependent hydrodynamic equations for a plane -parallel 
model of the solar atmosphere, using the method of characteristics. Radiative damp- 
ing by means of emission and absorption in H was taken into account. The numerical 
solutions gave the height where the waves steepen into shocks and dissipate most of 
their energy and the amount of energy that remains in the wave after shock formation. 

The temperature of two models of the solar atmosphere is shown as a function 
of the optical depth at 5000 A in Figure 1. The theoretical model in radiative equilibrium 
was constructed by Kurucz (1974). Its temperature decreases monotonically with height. 
The empirical model, the so-called HSRA (Gingerich etal. , 1971), was designed to 
reproduce the observations of the photosphere and the chromosphere. Its temperature 
drops through the photosphere to the temperature minimum and then rises again into 
the chromosphere and corona. The temperature at the minimum is approximately 
4000 K, where the optical depth is t 50 qq = 10 and the height above t 500 q = 1 is 
h = 550 km. 
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We have taken the static model atmosphere by Kurucz (1974) and have excited at 

its bottom a sound wave. The initial energy flux is that from Osterbrock 
7 2 

(3 X 10 ergs/cm sec), the wave period is 25 sec, the wavelength at the lower tem- 
peratures is 180 km, and the scale height is 120 km. The initial velocity amplitude 
is 0. 05 km/sec, and where the sound wave has steepened into a shock the wave velocity 
has reached 1 km/sec (Figure 2). Note that the sound velocity at that point is 7 km/sec. 
Thus, the shock is weak, and its Mach number is M = 1. 13. 

The shock occurs at a height of h = 650 km. This point defines the beginning 
chromospheric temperature rise and thus the location of the temperature minimum. 

The wave dissipates 90% of its energy at the shock and then continues with reduced 
amplitude, carrying the remaining 10% of its energy farther into the chromosphere. 

The temperature of the static atmosphere and the temperature perturbation caused 
by the wave scaled up by a factor of 10 are shown in Figure 3. At the shock front, the 
temperature perturbation reaches an amplitude of 500 K. 

The height of shock formation is shown as a function of the wave period and the 
initial energy in Figure 4. "HSRA observed" refers to the location of the temperature 
minimum in the empirical model. We see that, generally, the height where the shock 
forms increases with increasing period on wavelength, as expected. The increase in 
the height for short-period waves may be caused by the greater efficiency of radiative 
damping when the distance between wave maxima and minima is small; or we may be 
reading too much into these preliminary results. The "best" initial energy is 
3 X 10 ergs/cm sec. The temperature minimum for a 30-sec wave with this energy 
lies at 650 km, nearly one scale height higher than it does in the empirical HSRA. But 
this higher location may be in better accord with center-to-limb observations of solar 
radiation at long wavelengths. 

The wave energy carried beyond the point of shock formation is shown as a func- 
tion of wave period and initial energy in Figure 5. "HSRA" refers to the radiative flux 
emitted by the chromosphere in H - and in certain lines as determined by Athay (1970) 
and by others. We see that the agreement with observations is good for an initial 
energy flux of 3 X 10 ergs/cm sec and wave periods near 30 sec. 
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Figure 4. Height of shock formation as a function of wave period and initial energy flux. 
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Figure 5. Energy flux carried beyond the shock as a function of wave period and initial 
energy flux. 
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The main limitation of this work is that we have used only a single wave. Calcu- 
lation with several waves should give only minor changes in the main conclusions, 
however, since most of the energy is transported in waves with periods not very dif- 
ferent from the ones we used. 


CONCLUSION 

The low chromosphere is heated by short-period waves. In a calculation with a 

7 2 

single wave with a period of 30 sec and an initial energy flux of 3X10 ergs/cm sec, 

the wave dissipates most of the energy at a height (h - 650 km) that agrees with the 

observed location of the temperature minimum, and the energy flux (3 X 10 ergs/ 

2 

cm sec) that remains in the wave beyond the shock is sufficient to account for the 
radiative losses of the chromosphere. 
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FORMATION OF THE SOLAR EUV SPECTRUM 


Eugene H. Avrett 

This is a brief progress report on work in progress with J. E. Vernazza and 
R. Loeser on the determination of empirical chromospheric models. 

Figure 1 shows a theoretically determined continuum-intensity distribution at the 
center of the solar disk in the range 1300 to 100 A. This region is filled with emission 
lines, particularly in the wavelength regions longward of the indicated absorption edges. 
Only three emission lines are shown: the 304 A He II line, the 584 He I line, and Ly a. 
The insert shows the central part of the Ly a profile. 

In this brief report I will discuss only the continuum-intensity distribution, show- 
ing where in the atmosphere various spectral features originate, and how the spectrum 
is formed at those depths. The object of this work is to determine a temperature 
stratification, by use of a non-LTE model atmosphere in hydrostatic equilibrium, so 
that the computed spectrum agrees as well as possible with the observed one. After 
such a model is established for a particular type of solar region or structure, we can 
calculate the radiative and conductive contributions to the chromospheric energy 
balance. 

Before discussing the model, some comments should be made about how these 

o 

computed results agree with observations. The computed intensity at 1300 A and the 
Ly a wings agree reasonably well with the observed intensities, but the Ly a central 

o 

reversal is too deep. The computed C I ground state continuum at 1100 A is about a 
factor of 2 brighter than observed. The hydrogen Lyman continuum is in good agree- 
ment with the observations from 912 A to about 700 A. The He I 504 A continuum agrees 
approximately with quiet sun observations or is slightly higher. The calculated He H 
continuum is too bright, perhaps by a factor of 5 or more, but our He n calculations 
are very preliminary and have not yet been examined carefully. 
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Figure 1. 
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We are still in the process of adjusting the chromospheric model to match the 
observations, after having recently established a new model of the underlying photo- 
sphere and temperature minimum region (Vemazza, Avrett, and Loeser, 1975). 

Figure 2 shows the temperature distribution that was used to obtain these results. 
The temperature minimum of about 4200 K appears at the extreme right. The middle 
chromosphere is characterized by a broad temperature plateau at about 6500 K, 
between 1000 and 2000 km. These are heights above t = 1 at 5000 A in the photosphere. 
In the transition region between the chromosphere and corona we have introduced a 
20, 000 K temperature plateau in order to account for reversals in the higher Lyman 
lines and for the slope of the Lyman continuum. 

o 

Figure 3 shows details of the intensity calculation at 1300 A (or 0. 13 pm). Here, 

t is the monochromatic optical depth at this wavelength, dl/dh (on a linear scale not 

shown) is the contribution per unit height to the intensity I at p = 1, which at this 

wavelength is formed between the temperature minimum and the 6500° plateau. B is 

fib 

the Planck function, S is the continuum source function before scattering is taken into 

account, and S° is the source function after the scattering correction. J is the mean 

intensity. The scale on the left of the source function graph is the intensity in units of 
-2 -1 -1 -1 

ergs cm s st Hz The corresponding brightness temperatures are shown on 

the right. The upper left panel shows the relative contributions to the absorption 
coefficient while the lower left panel shows the relative contributions to the emission 
coefficient. The dominant processes are Si I, Cl, and absorption and scattering in the 
Ly a wing. The extreme departures from LTE are indicated in the right-hand panel by 
the differences between source function and Planck function. 

o 

Figure 4 shows the corresponding results at 1238 A in the far wing of Ly a. The 
place of formation is now at a greater height than before, and becomes much higher in 
the atmosphere as we approach line center. 


Figure 5 shows the results at the head of the 1100 A carbon continuum. As noted 
above this computed intensity is about a factor of 2 higher than observed. We can 
determine from this diagram where and by how much the temperature should be 
lowered to get a better result. 
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Figure 4. 











Figure 6 shows that the head of the Lyman continuum is formed at the other end 
of the 6500° plateau at the base of the transition region between chromosphere and 
corona. (The sharp variations in dl/dh and S C do not have any particular significance 
and should become smoother when more depth points are added. ) 

The higher temperature layers in the transition region become important at 

shorter wavelengths in the Lyman continuum, as shown in Figure 7. Here at 700 A 

~3 -1 

the contribution of the 20, 000° region between t = 10 and 10 is greater than the 

contribution closer to t = 1 . There is even a small contribution from yet higher 

-3 

temperature layers where t is less than 10 . This is an extreme example of the 

failure of the Eddington— Barbie r relation. 

Figure 8 shows that the conditions at the head of the He I 504 A continuum are 
strikingly different from those in the Lyman continuum. We have introduced an 
incident radiation field due to the various coronal emission lines, the strongest being 
the 304 A line of He n. The coronal radiation penetrates into the upper chromosphere 
causing He I ionization and excitation in a region where the temperature is 7000°. A 
local maximum occurs in the Dg and 10830 line opacities at around 2000 km. Reduced 
incident coronal radiation removes this extra component of absorption. As a con- 
sequence, these lines should be weaker, as they are observed to be, in the presence 
of coronal holes, as has been discussed in other papers at this meeting. 

Within the coming months we expect to complete a detailed study of the C I, H, 

He I, and He II continua and the H and He lines, and to establish the temperature and 
density stratification of quiet and active regions with different amounts of incident 
coronal line radiation. For each of these models we will determine the radiative and 
conductive contributions to the chromospheric energy balance. 
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